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HIGH  GAIN  FREE  ELECTRON  LASER  OSCILLATORS 

I.  Introduction 

We  have  conducted  an  analytical  and  numerical  analysis  of  the  field  evolution  in  a  hig-h 
gain  free  electron  laser  operating  in  the  oscillator  configuration,  as  depicted  in  Fig.  1.  The 
analysis  is  applicable  to  systems  with  electron  beam  pulse  lengths  which  are  longer  that  the 
particle  transit  time  in  the  resonator.  The  electron  beam  equilibrium  is  therefore  assumed 
to  be  spatially  uniform  and  temporally  stationary.  The  radiation  field  and  phase  averages 
which  are  performed  with  the  ensemble  of  electrons  is  conducted  for  an  interaction  length 
which  consists  of  the  entire  wiggler  structure.  This  is  in  contrast  to  other  simulations 
(theories)  which  perform  the  ensemble  average  over  the  wavelength  of  the  ponderomotive 
potential;  as  is  applicable  to  systems  with  temporally  stationary  fields1  or  short  beam 
pulses2  that  are  spatially  periodic. 

We  find  that  the  numerical  simulations  yield  qualitative  and  quantitative  agreement 
with  the  theory.  The  theory  for  the  example  given  (strong  pump  Compton  regime)  can 
be  separated  into  three  operation  regimes  which  we  shall  denote  as  the  ultra-high  gain, 
moderate  gain  and  low  gain  regimes.  Both  the  ultra-high  gain  (I\-L  >>  1)  and  the  low 
gain  (TkL  <<  1)  regimes  yield  growth  rates  that  exhibit  the  same  scaling  with  beam 
current,  energy  and  wiggler  field  as  is  obtained  for  an  FEL  amplifier  operating  in  these 
regimes.  Additionally,  we  consider  a  moderate  gain  regime  (T;X  >  1)  which  is  of  direct 
interest  to  NHL  experimental  parameters. 

II.  Theoretical  Model 

An  analysis  of  the  space  time  evolution  of  the  fields  and  particles  within  an  FEL  os¬ 
cillator  requires  a  self-consistent  coupling  of  the  fundamental  equations  for  the  particles 
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and  fields^5-8^  We  have  considered  a  Maxwell- V las ov  description  of  the  fields  and  par¬ 
ticles.  The  analysis  in  Appendix  A  results  in  the  following  system  of  equations  for  the 
coupling  of  the  fields  and  particles.  The  backward  travelling  wave  evolves  according  to, 
ddt{z)/dz  -  iadb(z)  =  0.  The  forward  travelling  potential  and  the  electrostatic  potential 
evolve  according  to, 


j-af(z)  +  iadf{z)  = 


-  Cl  f  dz'(z'  -  z)  exp[-iAA'(r'  -  z)\{d/dz'  -  iK)\3wdf{z')/2  -  o(z')\ 

Jo 

+  ic2  j  t/i'exp[-t'AA'(r'  -  .-)](a/3.-'  -  iA')((X  +  &)3wiI\z')l2  -  (1) 

Jo 


~o{z)  -  iKo(z)/ 2  = 


-  j  dz' -  r)exp[-iAAV  -  z)\(d/dz'  -  iK)[3.vdf(z')/2  -  of.-')] 


+  /  dz'  exp[-/AA'( -  z)](d/dz'  -  iK){Jwdf(z')/2  -  (1  -  J?0)o(*')J. 

J  o 


The  parameters  in  Eqs.  (1)  and  (2)  are  given  by,  a  =  Aw/c-  -:2(1  -  J*  )/2~0c*?.  ct  = 
c2(l  -  d?0)(~o  +  -M/"; o.  c2  =  ^3w/2~0c3;0~.  c3  =  ci^/3wR To.  c4  =  e3/(l  -  i?0), 
K  =  A,)  +  and  AA  =  A"  -  («^o  +  A~’)/i By  making  use  of  the  convolution  theorem, 
the  Laplace  transform  of  Eqs.  (1)  and  (2)  yields, 


{(.?  +  ia)[U  -  iAA')2+2c3]  +  iKci,3w/2}a/{s)  = 

{(3  —  iAA-)2  +  2cs(l  -  J2 /4)}d/(0),  (3) 

{(.9  -  f'AA')2  +  2c3}o[s)  =  L*.  -  id;(0)/A'|.  (4) 

where  d/($)  and  o(s)  are  the  Laplace  transformed  vector  and  scalar  potentials,  and  we  have 
retained  only  the  terms  in  the  driving  current  that  arise  from  the  momentum  derivatives  of 
the  phase,  DSKfdp;.  We  have  also  assumed  that  the  electron  beam  enters  the  resonator 
unbunched  so  that  o(r  =  0)  =  0. 

Since  the  singularities  of  the  Laplace  transformed  potentials  are  isolated  poles,  the 
Bromwich  inversion  of  these  transforms  can  be  easily  performed. 


<P{2)  =  2jResidue{0(s),sJ}exp(3J^),  (5) 

] 

afl2)  =  X^Residue  {«/(•»),  Sj}exp(S]3),  (6) 

i 

where  8j ,  are  the  poles  of  the  Laplace  transformed  potentials.  The  solution  for  the  radiation 
potential  given  in  Eq.  (6)  is  of  the  same  generic  form  as  the  solution  obtained  by  Bernstein 
and  Hirshfield9  for  the  FEL  amplifier  configuration.  Our  analysis  shall  differ  in  that  the 
backward  travelling  waves  are  not  neglected,  and  the  combination  of  the  forward  and 
backward  travelling  waves  are  required  to  satisfy  the  appropriate  boundary  conditions 
at  the  mirror  surfaces.  Specifically,  the  tangential  component  of  the  electric  field  must 
be  zero  at  the  non-transmitting  mirror  surface,  i.e.,  5/(0)  +a<,(0)  =  0.  At  the  partially 
transmitting  mirror  surface  at  the  end  of  the  resonator,  the  tangential  components  of 
the  electric  field  must  be  continuous,  i.e.,  5/(L)  exp(  —  ik^L)  +  5fe(I,)  exp(ikqL)  —  (1  — 
\/R)df(L)  exp(-iki)L)  =  0,  where  L  is  the  length  of  the  resonator  and  R  is  the  fractional 
power  reflected  from  the  far  end  resonator  mirror.  This  yields  the  following  expression  for 
the  boundary  conditions  at  the  mirror  surfaces, 


exp  — =  r^yexp{»(^3(l-j2/2)-2*o]L}^Residue{a/(s).s>}exp(a;L),  (7) 

which  is  the  equation  that  self-consistently  determines  the  complex  operating  frequency, 
A-;,  of  the  oscillator. 


III.  Results  for  Compton  Regime 

In  the  Compton  regime  the  effect  of  the  electrostatic  potential  can  be  neglected1 10,1 1\ 
In  addition  we  shall  assume  that  the  spatial  derivative  of  the  vector  potential  in  the  driving 
current  is  negligible.  These  derivatives  are  negligible  when.  j(s  -  *'AA')2j  >>  2 c3  and 
|j|  <<  K.  Under  these  conditions  the  Laplace  transform  of  the  vector  potential  is  given 


<*/(*)  =  5/(0)(s  -  1AA')2/  (a  -  Sj ), 


where  s}-  are  the  roots  of  the  dispersion  relation,  (a  -  «AA')2(a  +  /a)  =  -icq  3wK/2  ^ 
-*Wp(l  +  0zo)^lJkw/A^{)c2^.  Since  k0  is  a  free  parameter,  choose  k0  such  that,  AA'  = 

ko  3-  kw  —  (u^o  +  A«c)/tqo  =  -^/e  +  ^(l  -  ^/2)/2jJocrf.  Which  results  in  the  following 
solutions  to  the  dispersion  relation, 

A  .  .Ji  i2/y/z 

3j  =  ”*["T  "  2^(1  ~  ^/2)]  +  r»  1  -  ‘/n/3  .  (9) 

C  0  7  -1  -  i/y/Z 

where  T0  =  +  Jro)^^/4A^c27]1/,s/2Jz0  is  a  spatial  growth  rate  corresponding 

to  the  largest  spatial  growth  rate  in  the  amplifier  case12.  By  evaluating  the  residues  of 
d/(s)  for  each  of  these  poles,  one  obtains  the  following  solution  for  the  spatial  structure  of 
the  radiation  potential, 


*/(-)  =  5ZexP(s;J) 


It  is  evident  from  Eq.(lO)  that  the  spatial  growth  of  the  radiation  field  can  be  de¬ 
scribed  by  the  interference  of  three  modes;  which  can  be  identified  as  the  positive  and 
negative  energy  beam  inodes,  and  a  transverse  electromagnetic  inode.  The  constructive  or 
destructive  nature  of  this  interference  is  dependent  on  the  values  of  the  physical  parame¬ 
ters  which  characterize  the  roots,  .  For  physical  parameters  such  that.  I\)£  >>  1,  the 
unstable  mode  dominates  and  one  obtains  exponential  spatial  growth  at  the  rate  To. 

The  temporal  growth  rate  of  the  radiation  field  is  obtained  from  the’  negative  imaginary 
part  of  the  complex  oscillator  frequency,  A~.  The  oscillator  frequency  is  determined  by 
the  boundary  conditions  as  expressed  in  Eq.  (7).  which  for  the  approximate  roots  under 
consideration  vields. 


exp ~~~~  =  ~~?xp { ' [ ( 1  -  44)  -  2*. >11}  Vexp(,,L). 

c  O  *>  W  Qt  / 

J=1 
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We  shall  consider  three  distinct  solutions  to  this  equation  for  the  complex  oscillator  fre¬ 
quency.  The  first  of  which  is  the  ultra-high  gain  regime  (To£  >>  1).  in  which  case,  only 
the  fastest  growing  mode  in  Eq.  (11)  is  retained.  The  second  case  is  the  moderate  gain 
regime  (To-t  >  1),  where  only  the  decaying  mode  in  Eq.  (10)  is  neglected.  The  final  case 
is  valid  for  arbitrary  gain  and  all  terms  in  Eq.  (11)  are  retained.  The  imaginary  part  of 
the  oscillator  frequency  yields  the  following  temporal  growth  rates, 


1 ,  \/R  _  _  . 

-In  —  +r0I/2, 

\  In  ^  +  -  ln[2cosh(r0I)  +  2  cos(v/3r0I)], 

M  O  T 


Ultra-High  Gain 


Moderate  Gain 


iln  ^  +  iln[i  + 4008(^01) cosh(r0I) +  4 cosh2 (r0I)]. 

^  v  4 

Arbitrary  Gain 


In  each  of  the  expressions  for  the  growth  rate  the  first  term  is  negative  definite.  This 
represents  the  effect  of  losses  at  the  mirrors  and  che  coupling  losses  due  to  the  splitting  of 
the  radiation  into  three  modes.  The  necessary  condition  for  the  oscillator  to  lase  is  that 
the  remaining  terms  exceed  this  loss.  For  the  ultra-high  gain  case  this  requires  ToL  > 
-ln(\/^/3).  This  expression  has  been  confirmed  experimentally13  in  recent  operation 
of  the  XRL  FEL  oscillator.  The  interaction  length.  L,  can  be  varied  by  dumping  the 
beam  at  different  axial  locations  within  the  wiggler.  For  the  following  set  of  experimental 
parameters,  beam  energy  Eq  =  500  keV,  beam  current  I  =  100  .4,  wiggler  field  strength 
Bw  =  615  G,  beam  radius  r>,  =  0.64  cm  and  wiggler  length  l w  =  4.0  cm,  the  minimum 
interaction  length  is  determined  to  be  45  cm.  Inserting  this  value  into  Eq.(12)  yields  a 
theoretical  value  of  0.64  for  the  reflection  coefficient.  The  independently  measured  Bragg 
reflection  coefficient  has  the  value  0.65.  which  is  in  excellent  agreement  with  the  theoretical 


IV.  Multi-Mode  Simulation 


The  space-time  evolution  of  the  fields  in  the  resonator  is  simulated  by  numerically 
evolving  the  equations  for  the  fields  and  particles.  The  radiation  field  model  for  the  multi- 
mode  simulation  is  given  by, 

Afi(z,0  =  XX  ( t )  sin  (A:Bc)exp(i^’nf)el  +  c.c.,  (15) 

n 

o{ t)  =  y^Qin(t)  sin[(A*n  -f  kw)z  -  unt]  +  o2n(<)  cos[(An  +  kw)z  -  (16) 

n 

where  kn  =  nz/L  =  uin/c  and  the  sum  is  over  the  discrete  number  of  modes  under  con¬ 
sideration.  This  model  has  the  property  that  the  complex  expansion  coefficients  in  the 
harmonic  analysis,  an(f),Oin(f),0ora(f),  are  only  functions  of  time;  which  results  in  ordi¬ 
nary  differential  equations  for  the  particle  and  field  evolution.  This  model  also  has  the 
attribute  that  the  field  boundary  conditions  at  two  perfectly  reflecting  mirrors  is  automat¬ 
ically  satisfied,  and  we  model  the  resonator  losses  heuristically  by  adding  a  damping  term 
to  the  wave  equation,  \d2  jdz2  —  c~2d2/dt 2  —  i/c~2 d /dt]AR[z.  t)  —  izc~l  J±[z,  t).  where 
u  =  Q  and  Q  is  the  quality  factor  of  the  resonator.  The  driving  currents  and  charge  den¬ 
sities  for  the  vector  and  scalar  potentials  are  modeled  with  a  discrete  distribution  function 
as  follows, 


p(~,  t)  =  -e  J  dz0n,){t)fi{z  -  rU0.  01-  (17) 

J±{zJ)  =  [  dz0n0{t)6{:  -  z[z0.t))fr 0.  (18) 

me  J 

where.  A  =  A  R  +  Au;.  nn(0  =  «o(oc)[l  -  exp( -fjl-R )]  and  ri,.(  rc )  is  the  flattop  density 
of  the  electron  beam  pulse.  tR  is  the  characteristic  rise  time  for  the  beam  current  or  density, 
and  f(co-  0  is  the  axial  orbit  of  a  particle  located  at  position  at  t  —  0. 

The  slowly  varying  field  approximation,  dan[f)jdf  <<  *„an(/).  yields  the  following 
set  of  equations  for  the  evolution  of  the  fields  and  particles. 


(19) 


T» r  'r  .■  ■  'j-j  ■  1 !«  ■  .■  ■  ,■  ■  *  my  ■  y  ■  i  ■  w  ■  ?  ,  '.1  ■'/.■  ■  \  ■  »  ■  *  ■  i  ■  v  rHTT^ 


Aji  =  ~  3  in  /  dr0sin{yn(r0,r)]— P(r0), 

*  J  T  —  l  7o 

0'n  =  - [  dr0cos[e-n(r0.r)]— P(r0), 

*  An  ir-i  70 

= -.^n  /  dr0cos[ii-„(r0lr)]P(ro). 

Jt-1 

<4™  =  ~4n  /  c/r0  sin [tyn(r0,  r)jP(r0). 

J  T—l 

<  =  +  cos[(i>m(r0,  r)]  -  <?lm  sin[c’m(r0,  r)]] 

m 

~t”  *^4n  j  [(&m  "t*  sin[v^  (m)«  ^)] 

m  ^ 

-  cos|ftn(ro,r)j\. 


(20) 

(21) 


(22) 


(23) 


We  have  introduced  the  following  normalized  parameters,  r  =  Vqzt/L,  An  =  eAnjmc2 . 
O  =  eojmc 2.  (...)'  =  d[...)/dr.  We  have  also  made  the  following  definitions. 
dn  =  —  i An  (f  )exp(  i6n  r) ,  t-'rt('"0-t‘)  =  i^n  *b  '  )  •^jiF</*’0z  —  4(')i 

vL/v, 0z.  uj  =  4KLF[  1  -  d2/2)/2/:-c2T.  Ju  =  F  3wL3,r^ /2c2  320zkn.  3m  = 
2FJp30/c2{kn  +  M2  Jor,  3  in  =  L(kn  +  M /"7:4V  ^  =  4(*n  +  K)LF /232z-2 

and  P(r0)  =  1  -  exp( -tq/tr),  where  L  is  the  length  of  the  resonator.  F  is  the  filling  factor 
and  wp  is  the  nonrelativistic  plasma  frequency. 

This  system  of  equations  is  solved  numerically  by  using  a  four-point  Adams-Bashforth 
predictor  corrector  scheme  which  is  initialized  by  using  the  three  point  Runge-Kutta 
method.  The  ensemble  average  over  initial  electrons,  JT__l  (...)rfro.  is  typically  performed 
with  two  thousand  (2000)  particles.  The  results  of  the  simulations  and  the  linear  theory, 
obtained  from  the  linearization  of  Eqs.  (10)  -  (23),  are  shown  in  Figs.  2  and  3. 


V.  Conclusions 

A  comparison  of  the  temporal  growth  rates  obtained  from  the  linear  theory  and  the 
numerical  simulations  is  shown  in  Figs.  2  and  3.  The  growth  rates  for  the  simulations  are 
obtained  numerically  from  the  field  amplitude  data  during  the  initial  field  evolution,  where 

i 

I 


the  wave  growth  is  linear.  In  Fig.  2  the  data  is  presented  for  a  low  gain  case  with  physical 
parameters  given  by,  7  =  2.0,/  =  5.4.  F  =  0.2.  kwrb  =  0.62831.^*,  =  0.2  mdL/(w  =  50. 
There  is  an  excellent  agreement  between  theory'  and  simulation.  In  Fig.  2  we  also  compare 
the  theoretical  and  numerical  efficiencies  for  a  high  gain  case.  Where  the  efficiency  is 
defined  as  the  stored  electromagnetic  energy'  density  normalized  to  the  incident  beam 
energy  density.  The  theoretical  estimates  of  the  efficiency  are  based  on  particle  trapping 
arguments12,  with  the  assumption  that  all  the  energy'  lost  by  the  particles  is  converted 
into  electromagnetic  energy.  The  characteristic  change  in  velocity  of  a  particle  is  given  by 
the  difference  in  the  beam  velocity  and  the  phase  velocity  of  the  trapping  potential.  This 
phase  velocity  is  approximated  from  the  results  of  the  linear  dispersion  relation.  Again  we 
find  good  qualitative  and  quantitative  agreement  between  the  simulation  and  theory. 
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Appendix  A 


In  the  following,  we  shall  consider  the  space-time  evolution  of  the  radiation  fields 
produced  by  the  interaction  of  a  beam  of  relativistic  electrons  with  a  helical  wiggler  field 
contained  within  the  mirrors  of  an  optical  resonator.  The  analysis  is  fully  relativistic 
and  is  conducted  self-consistently  within  the  framework  of  the  Vlasov-Maxwell  system  of 
equations. 

The  wiggler  vector  potential  is  modeled  as  follows, 

Aw(s)  =  ^-[exp(«*w;)e_  +  exp(-i*wc)e+],  (.41) 

I+TJJ 

where  the  wiggler  magnetic  field  strength  is  Bw,  the  wiggler  period  is  iw  =  2*1  kw  and 
the  basis  vectors  are  e±  =  ( ex  ±  iey)j 2.  We  have  assumed  in  this  model  that  the  beam 
radius  is  small  compared  to  the  wiggler  period  (£*./•(,  <  1)  hence  the  transverse  gradients 
in  the  wiggler  field  are  neglected.  We  similarly  invoke  the  para-axial  approximation  to  the 
radiation  fields  and  neglect  transverse  coordinate  dependencies  in  the  fields,  to  obtain  the 
following  radiation  field  model  for  the  vector  and  scalar  potentials, 


A(z,  t )  =  [a/(;.  t)exp{-ik0z)  +  ab(z,  Oexp(</r0r)]exp(iu0f)e_  +  c.c.,  (.42) 

<?{z,  t)  =  d>(z,  l)exp[-i(k0  +  K)z  +  iu!0t}  +  c.c.,  (.43) 

where  «/(;,  t )  and  cib{z,t)  denote  the  forward  and  backward  components  of  the  wave  field 
respectively,  and  =  c^o  is  the  frequency.  These  field  coefficients  are  assumed  to  be 
slowly  varying  functions  of  space  and  time  compared  to  the  radiation  wavelength  and 
temporal  period.  The  slow  spatial  dependence  of  the  field  coefficients  is  expressed  by. 
I  Q~ldQ/d:  |<<  k,).  with  Q  =  o/(r,  l),  db(z,  t).o(z.  ()  and  the  slow  temporal  dependence 
of  the  the  field  coefficients  is  expressed  by,  j  Q~ldQjdt  ]<<  ^0. 

The  space-time  evolution  of  the  fields  is  governed  by  Maxwell's  equations,  which  can 
be  cast  in  the  form,  ( d'1  jdz 2  -  c~252/d<2).4(r,  t)  =  4 zc~l  J±{:,  t)  and  d2  jd:2o(:,t)  = 
-4 zp(z,t).  The  driving  current  and  charge  densities  are  obtained  from  the  appropriate 
moments  of  the  Vlasov  distribution  function.  The  Vlasov  distribution  function  is  evolved 


3 


according  to  the  equation,  {d/dt  +  [pz  j  m~;)d  j  d  z  -  e[E  +  (px  B)/m^c]  •  d/dp}g(z,  p.  t)  = 
0.  By  making  use  of  the  fact  that  the  canonical  transverse  momentum  is  an  invari¬ 
ant  of  the  motion  and  assuming  that  the  beam  is  cold  in  the  transverse  direction  (e.g., 
g(z,  Px,  Py,  pz,t )  =  g(z.pz,  £)<5(P2)<5(Pj,)  )  the  evolution  of  the  reduced  distribution  func¬ 
tion  is  governed  by, 


pz  d  ’  dd  e2 

m-fx  dz  [  dz  2m'! tC2 


where  me2 73*  =  {m2c4  -I-  t*2p2  +  e2(.4  •  .4)  }lt2.  Since  d(Aw  ■  A^.'j/dz  =  0,  the  equilibrium 
distribution  function  satisfies  the  equation,  {d /dt+(pz/m'i0)d  fd z}g^0^  {z.  pz.  t)  =  0,  where 
me270  =  {m2c4  +  c2p\  +  e2B2  /&2  For  long  electron  beam  pulses  we  shall  consider 
spatially  and  temporally  homogeneous  equilibria  given  by  g^(z,pz,  t )  =  3^(p’)-  To  first 
order  in  the  perturbed  fields,  the  evolution  of  the  linearized  distribution  function  is  given 


p-  d 


a  ,  r  r  ,1 


1  J7  +  — ~a~  =  ~  e~r  + - 7— (-4-  .4*) 

[df  m~iodz)  [  dz  m'i0c2  dz 

The  solution  to  the  linearized  Vlasov  equation  is  formally  given  by, 


rn'ioc2  dz1  vz  J  dpz 


Linearizing  the  wave  equations  for  .4  and  6  one  obtains. 


d 2  1  <9 2 


-4/T  f  -e2  - 


d\z.  t )  =  4zn0e 


By  making  use  of  the  slowly  varying  coefficient  approximations,  the  components  of 
Maxwell’s  equations  can  be  expressed  as: 

+  >-p-_’fwexp[-i(A^  -  woO]  f  dp2—g{l){:,pz;  t),  (A9) 

2w0C7  kw  J  7o 

|  j:  -  =  i-^Q-exp[i(A^  -  h?0f)l  J  dpzg{1){z,pz,t),  (.411) 

where  mc2~  =  {m2c4  +  c2p20  +  e2B2 /&2  J1/2,  a(n)  =  /  <fpr(7/7o)n3(0)  (ps), 
w2  =  4 jrn0e2/m,  j3w  =  eB^/m^c2^  and  AT  =  fc0  +  ku>-  Also  note  that  in  the  fre¬ 
quency  regime  for  resonant  interaction  of  the  radiation  field  and  the  beam  particles, 
(*,’o  ~  2  jz‘  zkwc)  the  driving  current  for  the  backward  wave  is  negligible.  We  shall  solve 
this  system  of  equations  in  the  time  asymptotic  limit  for  which  the  forward  wave  oscillates 
at  a  single  complex  frequency.  aj{zj)  =  5/(r)exp(/Aw/).  where  Aw  is  to  be  determined 
self-consistently  from  the  boundary  conditions  at  the  mirror  surfaces. 

The  ponderomotive  potential  in  terms  of  the  time  asymptotic  field  coefficients  is  given 


.4W  •  .4  =  — ^-«/(c)exp{— i[A';  -  (w0  -h  Aw)f]}  -f-  nonresonant  terms.  (.412) 

Retaining  only  the  phase  resonant  terms,  the  formal  solution  to  the  linearized  Vlasov 
equation  can  be  written  as  follows, 


jj(1)U,P*-,<)  =  f  dz'——^-exp{  -  i[A'  - 
Jo  P: 

{djdz'  .iK)[hg.a, o{-.')\. 


(wp  +  Aw), .  c)g[n1(p-1 
v*  J’  '  dpz 


(.413) 


Inserting  this  result  into  the  Vlasov-Maxweli  system  of  equations  one  obtains 
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^  -  -»} 


1  $w  /  i\  T/ 


WW-fJOl^M-p  (01  +  s£* 


/'ds'MC p{  -  ■AAV  -  -  ,A-)[|l-  ^)j,,a/(:,)/2  -  £(;')],  (.114) 

Jo 


JL  _  -  Jv  Vo(-)  =  (1  -  3?o)  (up  +A^’) 


2Kc2J  J?0 


f  d-V  -  «)exp{  -  <AA'(--'  -  *)}(d/3*'  -  iin[,4.5/(.-,)/2  -  o(~')] 

Jo 

^S=jr  J~  -  iAK(;'  -  *)}(»/».-'  -  >K)[ ,J„o/l  :')/2  -  (1  -  4,)®U')]- 


+ 


(-•115) 


where  we  have  defined  AA'  =  A'  -  (u-'o  +  The  previous  set  of  equations  yields  a 

dispersion  relation  which  we  shall  refer  to  as  the  complete  dispersion  relation.  A  simplified 
dispersion  relation  is  obtained  by  noting  that  in  the  momentum  integration  .  the  results 
are  most  sensitive  to  changes  in  the  exponent,  AA(r'  -  r).  Retaining  only  the  terms  in 
the  integration  by  parts  which  are  proportional  to  d±K/dp:,  one  obtains  the  following 
simplified  system  of  equations, 


■  A~' 


1, 


_  p  Jw  ( 1  ~  '^Jp  )  («*'o  A -o') 


2w0  n  d3 


*0 


d-{~  -  ;)pxp{  -  iAA’(;'  -  :)) 

c  Jo 


(-116) 


tA 


yp  (1  -  J?(1)  (•*' o  +  A-/') 


2AVT  ,J?o 
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•  ■  A'  ■  *  «.* 

.VvZ>-V-Y 


j  dz'{z' -  z)exp{  -  iAK(z’ -  z)}{— -  iK)[-£-af{z')  -  o(z')}.  (.417) 

In  both  the  cases  of  the  complete  and  simplified  set  of  equations,  the  equations  are  of  the 
convolution  type  and  can  be  solved  by  Laplace  transform  methods.  The  text  of  the  paper 
consists  of  a  detailed  analysis  of  the  simplified  set  of  equations  in  the  Compton  regime. 


